Superluminal Tachyonlike Excitations of Dirac Fermions in a Topological Insulator 
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We have considered a system of two topological insulators and have determined the properties of 
the surface states at the junction. Here we report that these states, under certain conditions exhibit 
superluminous (tachyonic) dispersion of the Dirac fermions. Although superluminal excitations are 
known to exist in optical systems, this is the first demonstration of possible tachyonic excitations in a 
purely electronic system. The first ever signature of tachyons could therefore be found experimentally 
in a topological insulator junction. 



Topological insulators (TIs), a new class of materials 
rich with new concepts and promises, have attracted con- 
siderable attention in the condensed matter physics com- 
munity (2| . In the bulk, the system is electrically insulat- 
ing, driven by the strong spin-orbit coupling present in 
the system. The three-dimensional (3D) TIs host a single 
Dirac cone in the surface states that was confirmed exper- 
imentally by angle-resolved photoemission spectroscopy. 
Dirac fermions are also present in graphene [3|. However, 
in the TIs, unlike in graphene, there exists only odd num- 
ber of non-degenerate Dirac cones with spin-momentum 
locking that results in helical Dirac fermions Q without 
spin degeneracy. This spin chirality of Dirac fermions 
prevents them from backscattering and localization 5]. 
That makes those systems ideal for spintronics applica- 
tions or for quantum computing [2]. Until now, most of 
the attentions have been heaped on the surface states 
of a single TI. Here we show that, the junction surface 
states of two TIs, in certain situations, exhibit superlu- 
minal (tachyonic) dispersion of Dirac fermions. Tachyons 
have eluded detection until now, despite diligent efforts 
by the particle physicists worldwide. However, as we 
have demonstrated below, it could perhaps be found in 
the present solid state system. 

We consider a junction between two topological in- 
sulators. The junction surface is described by z = 0, 
and we assume that the system of TIs is isotropic in 
X and y directions. Therefore, the surface states are 
characterized by the x and y components of the wave 
vector, k^ and fc , and the surface wave functions de- 
pend on the z-coordinate and decay in both directions, 
i.e., the positive and negative directions of the z axis 
[see inset in Fig. 1(b)]. We label the topological in- 
sulator at z < as TI-1, and the topological insula- 
tor at z > as TI-2. We have found that for a gen- 
eral variation of parameters of TI-2, the junction sur- 
face states exhibit one branch with unique tachyonic 
dispersion relation. Although, they are not yet found 
experimentally, these "Uberlichtgeschwindigkeitteilchen" 
(faster than light particles) discussed by Sommerfeld [G] 
in 1905, and many others [7:] since then, have always 



been vigorously pursued (e.g., in the case of the neutri- 
nos [si [g) b y the particle physics community for many 
decades |ig |. 

We assume that the electronic states of both TIs are 
described by the same type of low-energy effective 3D 
Hamiltonian [ll|, Il3 ; which has the 4x4 matrix form 
and can be expressed as 



^Ti = £(fc)- 



M{k)a^ — lA^a^dz 



^2^- "'a: 



^2^+'^x 



M{k)a^ +[A^a^d,_ 



(1) 
where a^ {i = x, y, z) are the Pauli matrices, dz = d/dz, 

k = (kx, ky) is a two-dimensional (2D) wave vector, k^ = 
k^ ± \ky, and 



e{k) = C^-D^dl + D^{kl 



ky), 



M{k) = M„ + B,dl^B,{kl + k. 



(2) 
(3) 



For a topological insulator of the type BijSeg, the four- 
component wave functions, ^, corresponding to the ma- 
trix Hamiltonian ([l} determine the amplitudes of the 
wave functions at the positions of Bi and Se atoms: 
(Bi*, Se^., Bii , Sci ), where the arrows indicate the direc- 
tion of the electron spin. In the case of BijScg TI, the 
constants in the Hamiltonian ([1]) are [i2|, A^ = 2.2 eV-A, 



^2 = 4.1 eV-A, B^ = 10 eV-A2, B^ = 56.6 eV-A^, 
Ci = -0.0068 eV, D^ = 1.3 eV-A^, D^ = 19.6 eV-A^, 
and Mo = 0.28 eV. 

The unique properties of the bulk Hamiltonian ^ is 
that, for a single TI, it can produce surface states with 
massless relativistic dispersion relation, E w v-pk, where 
Vp — A2\/l — {Di/BiY [13 is the Fermi velocity. In 
the case of two TIs, at the junction we expect a cou- 
pling between two surface states belonging to different 
TIs. Within a simple model which includes phenomeno- 
logical coupling between massless relativistic states of the 
two TIs, it was shown earlier that the properties of the 
junction surface states strongly depend on the relative 
sign of the Fermi velocities of the two TIs, i.e., on the 
relative sign of A2 for TI-1 and TI-2 15]. Here we show 
that for a realistic 3D model [Eq. ([T])] of the TI, one can 



observe new and unique features in the dispersion rela- 
tion of the junction states. 

In what follows, we study the junction surface states 
within the realistic 3D model of the TI. For two TIs, 
we assume that both TIs are described by the Hamilto- 
nian of the same type ^ but with different constants. 
To distinguish the constants corresponding to different 
TIs, we introduce superscripts (1) and (2) for TI-1 and 
TI-2, respectively. Following the general pro cedure of 
constructing the surface states of a TI [ij, [ll| , we first 
determine for each TI the general bulk solution of the 
Schrodinger equation of the form 'I' oc e"^ ^ e''^'', where 
m — 1 and 2 for TI-1 and TI-2, respectively. Substi- 
tuting this form of solution in the Schrodinger equa- 
tion, 'HTi{k,dz)'$ ~ E'i', we obtain a secular equation, 
det H^7^(fc,A(™))-S 



= 0, for each TI (m = 1,2). 
For each energy E, this equation defines four values 
of x'or\k,E), a = 1,...,4. Each x'^\k,E) is dou- 
bly degenerate, which finally generates eight wave func- 
tions for each TI, m = 1 and 2, *i™^(fc, £;) e^='"'^e''^'^ 
(see Ref. [11 [11), where s = 1,2, a = 1,...,4, and 
^s™ {k, E) is a four-component wave function. 

Our goal here is to determine the surface states local- 
ized at the junction of the two TIs. Therefore, out of 
the four values of A^™-' (for each TI), we choose only two 
values: for TI-1 (z < 0) with ReA^™^ > and for TI-2 

[z > 0) with ReAa < 0. After selection of these AL™ 
(q ~ 1,2), our junction surface states then take the form 

E^d'2*i!^e^='v*^^2<o 

s = l,2 
a=l,2 

This type of wave function determines the localized 
junction surface states. The energy of the junction 
state is found from the condition of continuity of the 
wave function, $(2:), and the corresponding current, 
[6Wj!l /5k^}<^{z), {k^ — idz) at the junction between 
the two TIs. Solution of this continuity equations de- 
termines the dispersion relation, E{k), of the junction 
surface states. 

We keep the parameters of TI-1 fixed as for Bi2Se3 (as 
given above) and vary the parameters of TI-2. Our re- 
sults indicate that for general variation of the parameters 
of TI-2, the junction surface states show one branch with 
unique dispersion relation that resembles the dispersion 
of the tachyons Jo\. In Fig. 1 we show the dispersion 
relation for the junction state when only one parameter, 
Ai, of TI-2 is different from TI-1. In this case there is 
only one type (tachyonic) of dispersion. The group ve- 
locity, Vg = dE{k)/dk, corresponding to the tachyonic 
dispersion becomes infinitely large at fc « fco [point B in 
Fig. 1(a)]. Similar superluminal dispersion was known to 
be present (theoretically) for propagation of the optical 
pulses through inverted two-level system of atoms [I6] 
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FIG. 1: (a) Dispersion relation of the junction states. The pa- 
rameters of the TI-2 are the same as the parameters of TI-1 
except for Ai — 2.2 eV-A. Near the point B the dispersion re- 
lation shows the tachyonic behavior with superluminal group 
velocity, (b) The locahzation length of the junction states 
shown in panel (a). Near the tachyonic point (point B) the 
junction states are strongly localized with localization length 
~ 10 A. At the edge of the tachyonic branch (near points 
A and C) the states become delocalized. Inset: Schematic 
illustration of the junction (red line) between two TIs. The 
coordinate system is introduced as follows: z-axis is perpen- 
dicular to the junction surface, and x and y are in the plane 
of the junction. 



and in metamaterial photonic crystals [I3] with folded 
bands. Superluminal propagation can also be observed 
for propagation of light pulses through a media with 
anomalous dispersion relation [18]. In all these cases, the 
tachyonic dispersion is achieved for propagation of light 
pulses through a specially designed medium. In contrast, 
our system consists of just electronic degrees of freedom 
and propagation of electronic excitations exhibits tachy- 
onic dispersion. If confirmed experimentally, this would 
be the first example of superluminal tachyonic dispersion 
in an actual semiconductor (electronic) system. 

In Fig. 1(b) we show the localization length in 
the z-direction of the junction states of the tachyonic 
branch. The localization length is defined as l^ — 

min ■^ l/|ReAs^ I f-. Clearly, the maximum localization 

of the surface states occur near the tachyonic point B 
with localization length « 10 A. These strongly local- 
ized states are less sensitive to the bulk disorder of the 
TIs, which can help with the experimental observation of 
the tachyonic states. Away from the tachyonic point the 
junction states become weakly localized and finally, near 
the ends of the tachyonic branch (points A and C) the 
junction states become delocalized. The actual distribu- 



surface states has the form 
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FIG. 2: The electron density along the z direction for the 
tachyonic branch shown in Fig. 1(a). The density is shown 
for the junction surface state at fc = 0.084 A~^. The arrows 
next to the lines are the directions of the electron spin. The 
junction between the TIs is at z = 0. 



tion of electronic density for the two (up and down) spin 
components, is shown in Fig. 2 for one of the junction 
states near point B [see Fig. 1(a)]. Both spin compo- 
nents are occupied with final non-zero spin-polazition of 
the junction tachyonic states. Exactly at the tachyonic 
point the junction state is spin unpolarized. 

The existence of a tachyonic branch in the dispersion 
relation of junction states means that the energy disper- 
sion becomes a non-analytic function of the 2D surface 
momentum, i.e., the group velocity v = h~^dE{k)/dk is 
infinitely large. Due to the analytical dependence of the 
Hamiltonian ([TJ on the wave vector k, the non-analyticity 
in the dispersion relation is possible only if the Hamilto- 
nian is non-hermitian. In our case the junction states are 
decaying states, i.e., $ ex e'^"^, where the real part of Ac 
is non-zero. For these states the Hamiltonian becomes 
non-hermitian and tachyonic branches are therefore al- 
lowed. The existence of tachyonic dispersion does not 
violate Einstein's causality principle. The reason is that 
the superluminal group velocity describes the propaga- 
tion of not a signal but an analytical wave packet. The 
propagation of a singularity, i.e., the signal, is not de- 
scribed by the group velocity, so there is no violation of 
Einstein's causality relation (see Ref. |16l[l8|). 

The tachyonic excitations can be described by an effec- 
tive 2D Hamiltonian. This effective Hamiltonian should 
have a 2 X 2 matrix form, which takes into account the 
electron spin degrees of freedom. In addition, the Hamil- 
tonian should also support the non-analytical tachyonic 
dispersion relation and should be non-hermitian. The 
tachyonic Hamiltonian can be constructed, for example, 
by introducing the imaginary proper mass in the Dirac 
equation |13|. In our case, the tachyonic Hamiltonian is 
achieved by introducing an imaginary Fermi velocity in 
the massive Dirac equation. More precisely, the effec- 
tive Hamiltonian which describes the tachyonic junction 
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ihv^k_ 



i/iWjfc^ 
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= \(^z 



ihvi{ak), (5) 



where Aq is the effective mass of the tachyons, and iwj is 
the imaginary Fermi velocity. This effective Hamiltonian 
produces a tachyonic branch with dispersion relation of 
the form 



E, 



Tach 



(fc)-± 



Ag 



h?v^k 



2 1,2 



(6) 



Therefore, k < kg = Aq/Hv^ and the group velocity at 
k ~ k^ becomes infinitely large. For the tachyonic branch 
shown in Fig. 1 the parameters of the effective Hamilto- 
nian (P are Ag = 0.313 eV and Wj = 5.7 x 10^ m/s. 

The wave function corresponding to Hamiltonian (5) 
with energy spectrum (6) has the following form 



^ 



Tach 



&'^/'^cosa{k) 
^ sina(fc) 



ic 



(7) 



where k 



\/k^ + k'^, (j) — cos ^ k^/k, and a{k) 



^o--ETaoh(fc) 

2An 



The corresponding direction. 



of electron spin is characterized by an angle 2a rela- 
tive to the z axis and is given by n = (n^,n ,n^) — 
(sin2asin(/), sin2acos(/), cos 2a). Therefore, for a given 
tachyonic state, the z-component of the electron spin is 
cos2a(fc), i.e., the state is spin polarized. At fc = fcg, i.e., 
at a singular point of the tachyonic branch, the angle 
a = 45'^ and the electron state is spin unpolarized. This 
behavior is consistent with the exact distribution of the 
electron density shown in Fig. 2, which illustrates finite 
spin polarization of the tachyonic state away from the 
singular point B (fc = fcp). 

The results shown in Fig. 1 illustrate the existence of 
tachyonic branches under variation of just one parameter, 
Aj^, of Tl-2. By varying the other parameters, we can 
introduce additional junction states with more than one 
tachyonic branch. If the signs of the constants A2 for two 
TIs, i.e., the signs of the Fermi velocities for two isolated 
TIs are the same, then usually we observe two tachyonic 
branches as shown in Fig. 3(a). If the signs of A2 are 
opposite then the junction states usually have one or two 
tachyonic branches and one massless relativistic Dirac 
branch [see Fig. 3(b)]. This massless relativistic branch 
was predicted earlier within a 2D model Hamiltonian of 
the junction states [15| . 

As the topological insulators with different Fermi ve- 
locities are already available in the laboratories, perhaps 
we could propose ways to detect the elusive tachyons 
in these systems. One possible physical manifestation 
of tachyonic dispersion relation in our system would be 
the display of a singularity in the time resolved measure- 
ments of 2D ballistic electron transport [20] along the 
junction between the TIs. Under the ballistic condition 
the electron wave packet propagates with a group veloc- 
ity, which is determined by the corresponding dispersion 
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FIG. 3: Energy dispersion curves for the junction states. The 
parameters of the TI-2 are the same as the parameters of TI-1 
except (a) A^ ' = 3.0 eV-Aand ylj = 5.1 eV-A; (b) Af ' = 3.0 
eV-Aand A-^ = —4.1 eV-A. Branches A and C correspond to 
the tachyonic dispersion, while the branch B describes mass- 
less relativistic Dirac fermions. 



relation. Within the effective model of tachyonic branch 
this group velocity is 



= n-^dETe.ch{k)/dk = ±- 



hvfk 



The actual singularity will be smeared due to a finite 
width in fc-space of the electron wave packet. 

In conclusion, we have shown here that for a general set 
of parameters for the topological insulator Hamiltonian, 
the surface states at the junction between two TIs have 
at least one unique tachyonic branch, which describes the 
propagation of Direc fermion excitations with superlumi- 
nal group velocity. Although the excitations propagating 
with superluminal velocity are known (theoretically) to 
exist in optical systems, we show here that such exci- 
tations can actually be realized in purely electronic sys- 
tems. In these systems, the tachyonic excitations are well 
localized at the junction between the TIs and are suscep- 
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FIG. 4: Schematic illustration of the ballistic propagation of 
tachyonic excitations (electrons with positive group velocity 
and holes with negative group velocity) through the junction 
layer between two TIs. At the point of diverging group veloc- 
ity, the ballistic travel time as a function of electron energy 
shows a cusp- like singularity (solid line). The finite width of 
the electron wave packet results in smearing of this singular- 
ity (dotted line). Point B corresponds to the singular point 
in Fig. 1. 



where the positive (negative) group velocity corresponds 
to electron (hole) excitation. Therefore the ballistic 
transport time through a finite distance, L^, is 



h = 



L^/v„ ex 



VW 



„ v.^ y .„y k^. At a point k = /cq, i.e., where the 
group velocity is infinitely large (point B in Fig. 1) the 
transport time becomes very small. At this point the 
transport time as a function of the energy would exhibit 
a cusp-like singularity. That means, at that energy the 
time of electron transport will have a sharp minimum. 



tible to direct experimental observation. It would indeed 
be a remarkable feat for condensed matter and the mate- 
rials sciences if the first ever signature of elusive tachyons 
is actually detected experimentally in a topological insu- 
lator junction. 
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[fl Electronic address: tapash@physics.umanitoba.ca 

[2] M.Z. Hasan and C.L. Kane, Rev. Mod. Phys. 82, 3045 
(2010); X.-L. Qi and S.-C. Zhang, ibid. 83, 1057 (2011). 

[3] D.S.L. Abergel, V. Apalkov, J. Berashevich, K. Ziegler, 
and T. Chakraborty, Adv. Phys. 59, 261 (2010). 

[4] D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, 
J. Osterwalder, L. Patthey, J.G. Checkelsky, N.P. Ong, 
A.V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, 
R. J. Cava, M.Z. Hasan, Nature 460, 1101 (2009). 



[5] P. Roushan, J. Seo, C.V. Parker, Y.S. Hor, D. Hsieh, D. 

Qian, A. Richardella, M.Z. Hasan, R.J. Cava, A. Yaz- 

dani. Nature 460, 1106 (2009). 
[6] A. Sommerfeld, Nachr. Ges. Wiss. Gottingen, pp. 201- 

235, February 25, 1905. 
[7] O.M.P. Bilaniuk, V.K. Deshpande, and E.C.G. Sudar- 

shan, Am. J. Phys. 30, 718 (1962); G. Feinberg, Phys. 

Rev. 159, 1089 (1967); E. Recami, J. Phys. Conf. Ser. 

196, 012020 (2009); O.M. Bilaniuk, tbid. 196, 012021 



(2009). 
[8] A. Chodos, A.I. Hauser, and V.A. Kostelecky, Phys. Lett. 

B 150, 431 (1985). 
[9] R. Ehrlidi, Am. J. Phys. 71, 1109 (2003). 
[10] E. Recami, Foundation of Physics 31, 1119 (2001). 
[11] C.-X. Liu, X.-L. Qi, H.J. Zhang, X. Dai, Z. Fang, and 

S.-C. Zhang, Phys. Rev. B 82, 045122 (2010). 
[12] H. Zhang, C.-X. Liu, X.-L. Qi, Xi Dai, Z. Fang, and S.-C. 

Zhang, Nat. Phys. 5, 438 (2009). 
[13] W.-Y. Shan, H.-Z. Lu, S.-Q. Shen, New J. Phys. 12 

043048 (2010). 
[14] B. Zhou, H.Z. Lu, R.L. Chu, S.Q. Shen, and Q. Niu, 

Phys. Rev. Lett. 101 246807 (2008). 



[15] R. Takahashi, S. Murakami, Phys. Rev. Lett. 107, 

166805 (2011). 
[16] R.Y Chiao, A.E. Kozhekin, G. Kurizki, Phys. Rev. Lett. 

77, 1254 (1996). 
[17] P.Y. Chen, C.G. Poulton, A.A. Asatryan, M.J. Steel, 

L.C. Botten, C. Martijn de Sterke, R.C. McPhedran, New 

J. Phys. 13, 053007 (2011). 
[18] P.W. Milonni, J. Phys. B: At. Mol. Opt. Phys. 35, R31 

(2002). 

[19] U.D. Jentschura, Arxiv preprint larXiv: 1201.63001 ('2012') . 
[20] E.A. Shaner and S.A. Lyon, Phys. Rev. Lett. 93, 037402 

(2004). 



